Soliton dynamics of an atomic spinor condensate on a Ring Lattice 



Indubala I. Satija 1 , Carlos L. Pando L. 2 and Eite Tiesinga 3 

1 School of Physics, Astronomy and Computational Sciences, George Mason University, Fairfax, VA 22030, USA 
2 IFUAP, Universidad Autonoma de Puebla, Apdo, Postal J-48, Puebla 72570, Mexico and 
3 Joint Quantum Institute, National Institute of Standards and Technology and the University of Maryland, 
100 Bureau Drive, Stop 8423, Gaithersburg, MD 20899-8423, USA 
(Dated: January 25, 2013) 

We study the dynamics of macroscopically-coherent matter waves of an ultra-cold atomic spin-1 or spinor 
condensate on a ring lattice of six sites and demonstrate a novel type of spatio-temporal internal Josephson 
effect. Using a discrete solitary mode of uncoupled spin components as an initial condition, the time evolution 
of this many-body system is found to be characterized by two dominant frequencies leading to quasiperiodic 
dynamics at various sites. The dynamics of spatially-averaged and spin-averaged degrees of freedom, however, 
is periodic enabling an unique identification of the two frequencies. By increasing the spin-dependent atom- 
atom interaction strength we observe a resonance state, where the ratio of the two frequencies is a characteristic 
integer multiple and the spin-and-spatial degrees of freedom oscillate in "unison". Crucially, this resonant state 
is found to signal the onset to chaotic dynamics characterized by a broad band spectrum. In a ferromagnetic 
spinor condensate with attractive spin-dependent interactions, the resonance is accompanied by a transition 
from oscillatory- to rotational-type dynamics as the time evolution of the relative phase of the matter wave of 
the individual spin projections changes from bounded to unbounded. 

PACS numbers: 03.75.Ss,03.75.Mn,42.50.Lc,73.43.Nq 



I. INTRODUCTION 

Spinor condensates are atomic Bose-Einstein condensates 
(BEC) with an internal spin degree of freedom that com- 
bine magnetism with condensation. Examples are optically- 
trapped spinor condensates of atomic rubidium-87|[l}Sl 5 
sodium||5]|6l, or chromium [7| with either spin or angular mo- 
mentum / = 1, 2, or 3 condensates with a three-, five-, or 
seven-component vector order parameter. In contrast to mix- 
tures of two or more atomic states] 8 1 or mixtures of several 
atomic species, spin-changing collisions in spinor conden- 
sates permit coherent dynamics among the hyperfine states. 
In a typical process for two / = 1 atoms, one in spin compo- 
nent m = — 1 and one in m = +1, can reversibly scatter into 
two atoms with spin component m = 0, which conserves the 
global magnetization of the condensate. This coherent spin 
mixing leads, nevertheless, to oscillations of the spin popu- 
lations, and is an analogue of Josephson oscillations in ultra- 
cold atoms l9l- tT3l . Hence, just as collisional interactions allow 
for a single-component condensate to be spatially coherent, 
spin-changing collisions, driven by internal interactions, allow 
coherence among internal degrees of freedom. A positive or 
negative sign of the strength of the spin-changing interaction 
determines whether the systems behave anti-ferromagnetic or 
ferromagnetic, respectively iflOlfTTl . Spinor physics has also 
been studied in optical lattices with exactly two atoms per lat- 
tice site fl5l . 

In view of the nonlinear nature due to interparticle inter- 
actions and high degree of control in the experiments, BEC 
systems are ideal systems for visualizing a wide variety of 
nonlinear phenomena. This includes solitary waves, the lo- 
calized nonlinear traveling waves that retain their shape, size 
and speed during propagation [16|. Experimental realization 
of solitons in a homogeneous single component BEC systems 
is a hallmark of the quantum coherence associated with many 



body systems. lfT71 - [T9l There has also been theoretical studies 
of solitons in a homogeneous spin-1 condensate [20 21]. Fur- 
thermore, theoretical and experimental studies of double-well 
bosonic Josephson junctions have unveiled novel phenomena 
such as broken symmetry macroscopic quantum self-trapping 
E2l and 7r-modes[23-25| as well as symmetry restored swap- 
ping modes[22|. In addition, numerical investigation of BEC 
systems, spatially separated into a ring lattice have demon- 
strated chaotic dynamics, deterministic dynamics with sensi- 
tive dependence on initial conditions. Il26l - l28l 

In this paper, we explore the quantum coherent time evo- 
lution of a spin-1 BEC that is spatially separated into six 
weakly-coupled sites arranged to form a ring geometry, see 
Fig. [T] For a large number of atoms such a system is well 
represented by a three-component wavefunction or order pa- 
rameter, $(a?, t), that satisfies a nonlinear Gross-Pitaevskii 
equation! 9, IQ)- Here, we restrict our calculations to the case 
with no external magnetic field. More importantly, we assume 
that all sites and spin components have at all times the same 
identical localized spatial-mode function, <f>(x). That is, for 
the m th -component of <£>(x, t) we have 

L 

^ m (x,t) = ^^(t)d>(x-x n ), (1) 

71=1 

where L is the number of sites, x n is the center of site n, and 
the dimensionless ^>™(i) are complex time-dependent ampli- 
tudes. The overlap between spatial-mode functions at differ- 
ent sites leads to tunneling between neighboring sites. 

More advanced numerical modeling uses time varying 
mode functions ||29l l30l or directly solves for the three- 
dimensional Gross-Pitaevskii equation lfT8l . As a first study 
in the soliton dynamics of a spinor condensate we believe that 
our simplified approach is justified. It shows the generaliza- 
tion of the spinor Josephson effect as well as the route to chaos 



2 




lattice index, 



FIG. 1 : (color online) Panel a) shows a schematic of a ring lattice of 
six sites (red dashed circle with labeled markers) containing a / = 
1 spinor condensate. The solid black line represents the superfluid 
population along the ring as a polar graph. Panel b) shows the initial 
fractional population, p n , as a function of site index n used in our 
dynamical simulations. The distribution corresponds to an excited 
(period-6) solution of the DNLSE. For each site the population of 
the three magnetic sublevels of the spinor is the same. 



in its clearest form. It should be noted that the ring lattice with 
less than three sites is an integrable system. BEC on a lattice 
with three or more than three sites exhibits amazing degree of 
complexity including coexistence of regular and irregular be- 
havior that is known to accompany chaotic dynamics. In our 
detailed numerical exploration of ring lattices with different 
numbers of sites, the key characteristics of the time evolution 
as discussed in this paper are independent of the number of 
sites. Here, we will present our results for a ring lattice of six 
sites only. 

Our initial state is a discrete solitary wave that is an ex- 
cited eigenstate of the discrete nonlinear Schrodinger equa- 
tion (DNLSE) satisfied by the ip™ (t) in the absence of spin- 
changing interactions. The spin-changing interaction then 
induces oscillatory dynamics, an internal Josephson effect 
where different spin components evolve with different ampli- 
tudes and phases. In reality the spin-changing Hamiltonian 
can not be turned off and this state needs to be engineered 
by a combination of resonant electro-magnetic radiation that 
controllably creates superpositions of m-states [6] and off- 
resonant spatially-dependent light forces that induce solitons 

nana. 

In this article we show that for weak spin-changing interac- 
tions the dynamics on the ring lattice is dominated by two 
frequencies. Interestingly, one of these frequencies, which 
we denote as ui, describes the periodic behavior of spatially- 
averaged spin populations. The second frequency, denoted by 
17, describes the dynamics of the site-dependent population 
summed over the three spin components. In other words, the 
dynamics of collective or averaged coordinates is pendulum- 
like. As we discuss below, this is somewhat surprising as the 
spatial profiles of the three spin components evolve differ- 
ently and the dynamics on the ring should deviate from that 
for a "simple" trapped spinor condensate. The single-mode 
approximation for the latter system has been shown to lead 



to pendulum-type physics iflOl [T2ll . Our study suggests that 
the origin of this "regular" global dynamics has its roots in 
the strong correlations among different lattice sites and spin 
projections and is not rooted in thermal averaging over many 
degrees of freedom. We note that u is solely determined by 
the spin-changing interaction, whereas £1 shows a rather weak 
dependence. We find that f2 is mainly controlled by the size 
of the ring. 

By increasing the absolute value of the spin-dependent in- 
teraction strength, the two frequencies can be mode locked. 
This resonance condition is found to describe the onset to 
chaotic dynamics. In other words, the spinor condensate ex- 
hibits a transition from quasi-periodic dynamics to chaotic 
dynamics where at the onset to transition, all local and 
global degrees of freedom oscillate in unison. Consequently, 
the initial soliton profile reappears periodically, providing a 
unique demonstration of a novel-type of spatio-temporal in- 
ternal Josephson effect. Furthermore, in contrast to the anti- 
ferromagnetic case, in the ferromagnetic condensate, the res- 
onance is accompanied by a transition from bounded to un- 
bounded dynamics. In the anti-ferromagnetic case, the dy- 
namics remains bounded and oscillatory. Nevertheless, be- 
yond the resonance additional frequencies appear. 

For sake of simplicity, we restrict our simulations to zero 
magnetization. This is enforced by choosing initial condi- 
tions where the m = +1 and m = — 1 components have 
same initial wave function. In fact, we choose an initial state 
where all three components have the same initial wave func- 
tion. An alternate initial state with zero magnetization and 
contains mostly m = ±1 atoms lead to similar results. It will 
not be described here. 

In section [TT] we review the mean-field equations that de- 
scribe the evolution of the spinor condensate on a lattice as 
well as the underlying single-mode approximation approxi- 
mation on each lattice site. The latter has been shown to pro- 
vide a reasonable description for the continuum system. For 
our ring lattice, the SMA suggests collective coordinates to 
describe the global dynamics. The initial state is described 
in Sec. [Ill] In section |IV] we show numerical simulations 
for weak positive, anti-ferromagnetic spin-dependent integra- 
tions and show that the dynamics, although complex, is dom- 
inated by two frequencies. In sectionJV] we discuss the reso- 
nance condition that occurs when the spin-dependent interac- 
tion strength is increased. At resonance we observe an onset 
to chaos. Section [Vl|brie fly illustrates the ferromagnetic case. 
We conclude in Sec. I VIII 



II. MEAN-FIELD EQUATIONS FOR A SPINOR 
CONDENSATE AND COLLECTIVE COORDINATES 

An atomic spinor condensate with large atom number is de- 
scribed within mean-field theory by a complex vector order 
parameter or condensate wave function \&(x, t), whose evolu- 
tion is governed by a multi-component Gross-Pitaevskii equa- 
tion J9] [Toll . We will assume that the order parameter can 
be approximated by Eq. [JJ where a single time-independent 
mode function, <fr(x), determines the spatial dependence in 
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each well. We denote this by the L-site single-mode approx- 
imation (L-SMA) in analogy to the single-mode approxima- 
tion (SMA) for a spinor condensate in a dipole trap. 

The interactions between two spin-1 atoms have 
a spin-independent and a spin-dependent or (spin- 
changing) contribution. Within a mean-field the- 
ory and the L-SMA their strength is given by 
c q = Airh 2 g q /{2p) f dx\(f>(x)\ 4 / fdx\(j)(x)\ 2 with q = and 
2 for the spin-independent and spin-dependent contribution, 
respectively. Here p is the reduced mass for two atoms 
and h is the reduced Planck constant. The lengths g q are 
.9o = ( a o + 2a2)/3 and g% = (a 2 — ao)/3, where clq and a 2 are 
scattering lengths for s-wave collisions of two / = 1 bosons 
with total angular momentum F = and 2, respectively. 
For stable condensed gases we require that c > 0. We 
note that c 2 > for anti-ferromagnetic Na and c 2 < for 
ferromagnetic 87 Rb. 

The dynamics of the ijj™ (i) are governed by the DNLSE 

m^- 1 = -j(^-h 1 + ^-l 1 ) (2) 
+( Co + C2 )(|Va; 1 i 2 + i^i 2 )^ 1 
+( Co - C2 )l^i 2 c 1 + C2 (^) 2 (^r 
rn-^l = - J«_! + i>° n+1 ) + c |v£| V° 0) 
+( Co + C2 )(|^| 2 + icT)^ + 2c 2 ^v I ; 1 (^)* 
= - J(v£_i + ^ +1 ) + (c + c 2 )(|V^I 2 + l^l 2 )^ 

+ (c - c 2 )^" 1 | 2 ^ + ^(OW)* , (4) 

where J is the positive site-to-site tunneling energy and we 
use periodic boundary conditions. In the absence of the last 
term on the right hand side of Eqs.|2]|4] the set of equations de- 
scribes a three-species condensate. Spin-changing terms make 
a spinor condensate unique as they induce population oscilla- 
tions between m levels. Equations 2|4 conserve total atom 
number and magnetization. In other words, J2nm Wn W| 2 
and J2nm TO I' ! /'™(^)| 2 are conserved. 

We will monitor the local population of each spin state as 
well as global or collective coordinates, such as spatially- and 
spin-averaged population and phases. It is therefore conve- 
nient to define ip™(t) — *J p™(t) exp[i0" l (i)] with popula- 
tions p™(t) and phases 0™(t). Following Refs. J9j [TO | for a 
single-mode spinor condensate natural local canonical coordi- 
nates are 

Zn = P° n -(p+ 1 +Pn 1 ) and 7 n = 0- 1 +^-2$. (5) 

as well as globally-averaged coordinates Z = J2 n Z n /L and 
7 = *}2 n ^ln/L. Throughout this article we call 7„ and 7 
spinor phases. Other useful population averages are 



and ^ = 3 E Pn 

n—1 m— — 1 



(6) 



In a dipole trap, or equivalently for L = 1, a simple spinor 
model is given by the SMA. The canonical variables Z and 7 
then satisfy the pair of equations 



which are independent of the spin-independent interaction 
with strength cq. For small Z and 7 the dynamics are har- 
monic. In general, however, these coupled equations describe 
a nonlinear pendulum whose length depends upon the mo- 
mentum. 

For a L > 2 ring of lattice sites the collective variables Z 
and 7 satisfy 

h 



(8) 
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" 2 f zZ z ^ 1+cos ^ + 'i w 



(9) 



where W = MEnM+Jfn-^li/^-^nli/^n 1 )}- 
The spin-independent interaction strength Co does not explic- 
itly appear in these equations. As already mentioned in the 
introduction and further discussed in next section, this sug- 
gests that the spatially-averaged Z and 7 will exhibit periodic 
oscillations, with frequency oj, that are solely determined by 
the spin-changing interaction. Finally, the spin-averaged pop- 
ulations at each site satisfy 



J 



Im £ (Olff+i + C-ll . (10) 



£2 
H 



f'2 



(l-Z 2 )smj and 7 = -2-fZ(l +COS7) , (7) 



which does not explicitly depend on c 2 . Hence, we expect 
that the a n oscillate periodically, characterized by frequency 
Q, and its harmonics. 



III. INITIAL DISCRETE SOLITON 

For a ring lattice of six sites, we use the initial condition 
shown in Fig. [T] Commonly referred to as a discrete soliton, 
it is a stationary solution of the DNLSE in the absence of the 
spin-changing interaction [26 , 3T) (i.e. with c 2 = 0) and Co = 
J. Along with phase </>™ = or tt for even or odd site index 
n, respectively, the initial state populations are the same for 
the three components. 

This stationary solution, corresponding to a solution 
V'JTC*) — iftnit ~ 0) exp(— ipt/H) where p is the chemi- 
cal potential, is obtained by numerically solving the result- 
ing nonlinear map, as Eqs. |2|4] reduce to a set of L two- 
dimensional cubic maps. Following Ref. [28], we find that for 
L = 6 the solutions are 6-fold degenerate and for a range of p 
the localized soliton mode in Fig.[T|is the only stable solution. 
In fact, we have used p = 2.5 J. Intriguingly, for these values 
of the chemical potential the homogeneous solution with the 
same density at all sites is unstable. 

It should also be noted that known localized soliton-type 
solutions with zero phase for all sites correspond to attractive 
spin-independent interactions with cq < 0. To obtain local- 
ized solutions for repulsive interactions, one needs to consider 
solutions with phases different from zero. In general, it can be 
shown that for a lattice with even number of sites, the mapping 
ip™ — > (— l) n ('0™)* relates a soliton solution for attractive in- 
teractions with those with repulsive interactions. For a single- 
component condensate discrete solitons have been studied ex- 
tensively for ring lattices of various sizes 



FIG. 2: (color online) Time dynamics of an anti-ferromagnetic spinor 
soliton on a six-site ring for several spin- and/or spatially-averaged 
degrees of freedom assuming a small positive spin-changing inter- 
action energy. Time is in units of hj J, where J is the tunneling 
energy between the sites. Calculations are performed for C2 = 0.1J 
and Co = J ■ Panel a) shows the dynamics of the spatially-averaged 
spinor population Z (black line) and the spinor population Z„ for 
the individual sites n (colored lines). The symmetry of the initial 
soliton around n — 3 implies that only four of the six sites have a 
distinct time evolution. Panel b) shows the spinor phases 7 (black 
line) and "/„ (colored lines). Finally, panel c) shows the site-specific 
population a n , averaged over the three spin components. 



FIG. 3: (color online) Power spectrum or Fourier transform of 
the time evolution of the spatially-averaged spin population p m =+i 
(red), spatially-averaged spinor phase 7 (green), and spin-averaged 
population a n =i at site n = 1 (blue). The parameters and initial 
state are as in Fig. [2] and the frequency is in units of J/h. The top 
and bottom panel show the power spectra on a linear and logarith- 
mic scale, respectively. Also indicated are the dominant frequencies 
ui and Q, which have been assigned as due to the spin-dependent 
and spin-independent interactions, respectively. In the bottom panel 
higher harmonics of Q can be observed. 



Experiments with single-component Bose condensates in 
double-well potentials |24| have observed Josephson oscilla- 
tions and quantum self-trapping in the limit J < Co. The 
opposite limit can also be reached leading to tunneling of 
(nearly-)independent atoms. Here, we chose a compromise 
with Co = J. As an aside we note that with our initial state 
and c = J we have implicitly specified the chemical poten- 
tial and thus atom number in each site. 



IV. QUASIPERIODIC DYNAMICS OF AN 
ANTI-FERROMAGNETIC SPINOR 

Figure|2]shows the time evolution of the spinor soliton on a 
six-site ring, described in terms of the spinor coordinates Z n 
and Z, spinor phases j n and 7, as well as populations a n . We 
use a small positive spin-changing interaction strength c 2 for 
an anti-ferromagnetic spinor. We observe that the time depen- 
dence of Z and 7 are nearly sinusoidal. The site-dependent 
Z n , 7„ and a n , however, oscillate at a higher frequency. They 
do so in a non-sinusoidal manner with sharper minima than 
maxima or vice versa. The phases j n only show small excur- 



sions around the average 7. Finally, we note that the spinor 
phases are bounded for oscillatory motion. 

Figure [3] shows the power spectrum of three of the time 
traces shown in Fig. [2] It highlights the existence of two 
dominant frequencies, w and 57, as well as weaker higher 
harmonics in il indicating non-sinusoidal periodic behavior. 
The spatially-averaged degrees of freedom predominantly os- 
cillate with a frequency to, which from simulations with other 
small C2 is found to be proportional to the absolute value of c 2 . 
In contrast, spin-averaged local populations oscillate nearly- 
sinusoidal with frequency ft, which from other simulations is 
found to weakly depend on the spin-changing interaction but 
is inversely proportional to the number of lattice sites. Fi- 
nally, the local dynamics for individual spin projections is 
quasiperiodic with a slow frequency u> and a faster beat fre- 
quency at multiples of fl This illustrates the correlations that 
exist among the spatial and spin components. Detailed stud- 
ies with various ring sizes indicate that the spinor dynamics 
is characterized by two frequencies, irrespective of the size of 
the ring. 
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FIG. 4: (color online) Time series for an anti-ferromagnetic 
spinor showing the resonance condition for increasing spin-changing 
strength C2. Time is in units of h/ J and Co = J. In all panels Z (red 
curve) and <r n= i (blue curve) are shown. Panels a), b), and c) show 
results for C2 = 0.5J, 0.65J, and 1.0J, respectively. The resonance 
condition occurs for C2 = 0.65J, where the spatially-averaged Z 
and spin-averaged a n =\ oscillate at the same rate. For C2 > 0.65J 
the traces are chaotic. 



V. RESONANCE CONDITION FOR THE 
ANTI-FERROMAGNETIC CONDENSATE 

Figures [T|and|2]showed that for small C2 /cq two frequencies 
u and fi are very distinct. For increasing positive C2 at fixed 
Co, as shown by time traces in Fig. [4] and power spectra in 
Fig. |5J both frequencies increase although at a different rate. 
For a critical spin-changing strength C2 when uj = ft the anti- 
ferromagnetic spinor reaches a "resonance" state. All three 
spin-components at all the sites of the ring then oscillate as a 
single entity. For cq = J this resonance occurs at = 0.65 J. 
For large values of C2 the behavior becomes chaotic, which is 
apparent as broad-band features in the corresponding power 
spectrum shown in Fig. [5] 

Fourier analysis of collective as well as local coordinates 
shows that the resonance state corresponds to matching of not 
only the dominant frequencies cj and f2, but also some of the 
other less prominent frequencies (not visible in the linear plot) 
as illustrated in the log-linear plot on the middle row of Fig. [5] 
In fact, at resonance, the dominant or the primary peak is ac- 
companied by secondary satellite peaks, equally spaced on the 
either side of the primary peak. 

Figure [6] further illustrates that at resonance all sites oscil- 
late in phase with same frequency. The spatial profile of the 
soliton reemerges periodically without any significant change 
from the initial profile. Thus, the resonance state in an ordered 
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FIG. 5: (color online) Power spectra showing the resonance condi- 
tion for increasing spin-changing strength C2 for the same parame- 
ters as in Fig. [4] The frequency is in units of J/h. In all panels the 
spatially-averaged p m =+i (red curve) and spin-averaged a n =i (blue 
curve) are shown. Panels on the left and right show spectra on a 
linear and logarithmic scale, respectively, while from top to bottom 
c 2 = 0.5J, 0.65J and 1.0J. 




FIG. 6: (color online) Spin-averaged soliton profiles, <r n , for the anti- 
ferromagnetic spinor before and at resonance for five times spaced 
by time interval 2n/Q, starting with the initial state (red curve). We 
use Co — J and panels a) and b) correspond to C2 = 0.5J and 
C2 = 0.65J, respectively. The soliton reforms close to its initial 
state at resonance when C2 = 0.65 J. 



state where wave functions of all the components of the spinor 
condensate at all L-sites of the ring lattice oscillate in unison 
and the dynamics is well characterized one frequency. The 
fact that the soliton profile reemerges periodically provides a 
unique demonstration of quantum coherence and a novel-type 
of internal spatio-temporal Josephson effect. 



6 




300 



VI. FERROMAGNETIC CONDENSATE 

We now briefly describe the dynamics of a ferromag- 
netic spinor condensate with negative C2. For a small spin- 
changing interaction, the dynamics is similar to that of an 
anti-ferromagnetic spinor. Namely, the spatially-averaged be- 
havior is pendulum-like with dominant frequency u>, while the 
local population oscillates with frequency Q. An example is 
shown in Fig. [7] The local density for each spin component 
oscillates with both frequencies. A resonance state can again 
be achieved by increasing \c2\ but, for cq — J, now occurs 
when uj/fl = 3. 

In contrast to the anti-ferromagnetic case, however, the res- 
onant transition is accompanied by unbounded dynamics or 
phase winding where the phase of the condensate becomes un- 
bounded as shown in Fig. [8] This behavior manifests itself as 
a zero-frequency mode in a power spectrum. In other words, 
chaotic dynamics with its broad-band spectrum is accompa- 
nied by a transition from an oscillatory to rotational mode for 
the collective degrees of freedom. 



FIG. 7: (color online) Time evolution of spatially-averaged popula- 
tion Z (red) and spin-averaged population a n= \ (blue) for a ferro- 
magnetic spinor condensate with negative C2. Panel a), b), and c) 
show traces for C2 = — 0.075J, — 0.09J, and — 0.095J, respec- 
tively, corresponding to cases before, at, and beyond resonance. We 
use Co — J and time is in units of h/ J. 
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FIG. 8: Phase portraits or parametric plots of the time-evolution of 
the spinor population Z and phase 7 for a ferromagnetic spinor at 
C2 = — 0.075J, panel a), and C2 = —0.095 J, panel b). The system 
is at resonance in panel b). The dynamics in panel b) is accompanied 
by unbounded motion of the spinor phase. We use Co = J. 



VII. DISCUSSION 

In summary, numerical explorations of mean-field equa- 
tions of spinor condensate on a ring lattice with a small 
number of sites reveals strikingly correlated dynamics of the 
many-body system. Even though the condensate is separated 
into L-sites, soliton dynamics at all sites can be characterized 
by just two frequencies. With quasiperiodic dynamics at lo- 
cal sites for individual spin components, spatially averaged 
behavior for each spin component as well as spin-averaged 
dynamics at each sites is found to be periodic. The fact that 
the time series describing local dynamics is characterized by 
two frequencies and these two frequencies untangle in the col- 
lective degrees of freedom is rooted in correlations among dif- 
ferent spin degrees of freedom at various sites of the lattice. 
However, proper understanding of these correlations remains 
an open challenge. Our study with ring lattices of various 
sizes show that two-frequency characterizations of the spin-1 
condensate on a ring lattice is valid irrespective of the number 
of sites on the lattice. 

Our study provides a new illustration of a quasiperiodic 
route to chaotic dynamics in a many-body system where the 
critical point is known to be characterized by a resonant 
state. Simple models of dynamical systems such as the one- 
dimensional circle map ll32l are paradigms of quasiperiodic 
route to chaos, where the critical point corresponds to pa- 
rameters where the two frequencies are mode-locked. The 
quasiperiodic route to chaos is a well-established scenario 
in dynamical systems exhibiting a transition from regular to 
chaotic dynamics [32 1. However, the fact that the critical point 
describes periodic dynamics is a novelty in many-body sys- 
tems rooted in the coherence associated with spinor conden- 
sate. In this case, the critical point is a highly-ordered many- 
body state exhibiting a new type of spatial-temporal coher- 
ence. At the critical point where the two dominant frequen- 
cies are in resonance, the ring lattice with L-sites oscillates 
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in unison with a single characteristic frequency. This internal temporal coherence. 

Josephson effect where an unscathed soliton profile reemerges This research is supported by Office of Naval Research, the 
periodically provides a novel illustration of both spatial and CONACYT-Mexico and the US Army Research Office. 
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